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ABSTRACT 



Aims. The first (/-band and new BV(RI)c photometry of the RRd star V372 Ser is presented to determine some parameters of the star. 
Methods. In April, May 2007 2812 U, B, V,Rq, Ic frames were obtained at Konkoly and Teide Observatories, 1508 V observations 
were collected from the literature. Fourier fitted light curves have been derived in all bands. 

Results. The non-linearly coupled frequencies /„ = (2.121840 ± .000001) cycles • day -1 , f = (2.851188 ± .000001) c • d -1 , i.e. 
periods P = 0.4712891 ± .0000002 days, P t = 0.3507310 ± .0000001 d, Pi/P = 0.7441950, amplitudes = 0.15399 mag, 
Aj = 0.20591 mag, and phases have been found. Ai/Aq = 1.319 ± .008 has been found from averaging the amplitude ratio in the 
different bands i.e. the first overtone is the dominant pulsation mode. From the V observations upper limits are given for secular 
change of the Fourier parameters. The period ratio and period put V372 Ser among the RRd stars of the globular clusters M3 and IC 
4499, mass, luminosity, and metallicity estimates are given. 
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1. Introduction 



stars: fundamental parameters: (colours, frequencies, mass, luminosity) - stars: individual: 



Double mode pulsation plays a key role in the mass determina- 
tion of RR Lyrae (RRL) variable stars because confirmed RRL 
type component is not known in binar y syste ms. We have indi- 
rect methods to determine their mass (Smith]| 1995l): for d ouble 
mode RRL (RRd) stars the Petersen diagram dPetersenl 19731) of- 
fers a possibility of most accurate determination. Of course, the 
derived values depend on theoretical input. 

In spite of the importance of RRd stars there are very few 
observed data of their multicolour photometric behaviour. It is 
partially due to the moderate number of known stars until the 
near past: at the moment of this writing we know 30 more or less 
bright RRd stars (V < 14 mag) in the Galactic field discovered 
by th e NSVS and ASAS all-sky va r iability surveys dWils et alj 
120061: ISzczvgiel & Fabrvckvl 120071: iPilecki & Szczvgiell 120071 
lKhruslovl20d7h . Only two stars (GSC 4868-0831 and V372 Ser) 
are known with V < 11.5 mag. A n umber of fainte r objects 
were discovered in the Galactic bulge (Mizerski 2003) and sev- 
eral more in extragalactic systems allowing to draw statistically 
significant conclusions concernin g distribution of the RRd stars 
in the Petersen diagram, etc (see IClementini et"aT1l2004l and ref- 
erences therein.) 

To review the literature we found that AQ Leo, CU Com, 
two fainter Galactic, an d nine extragalactic RRd stars were 
studied spectroscopically dClement et al.ll99lt IClementini et al.l 
2000; iGratton et al]|2004l) . Multicolour photometric time serie s 
are available for CU Com (CCD BVI, IClementini et all [2000). 



Send offprint requests to: S. Barcza 

* Photometric data are only available in electronic form at CDS 
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BS Co m (CCD B VRI, lDekanvl|2008l) and AQ Leo (photoelec 



nic BV, Jerzvkiewicz & Wenzelll977tlJerzvkiewicz et aljl982b . 
Some multicolour CCD time series were also published for six 
globular clusters containing at least one RRd star (IC 4499, M15, 
M68, M3, NGC 6426, and NGC 2419). 

The RRd character of V372 Ser (GSC 5002-0629, a = 
15 h 17 m 35 5 .03, 6 = -l°5'17'/ 2) was discovered a few years ago 
dGarcia-Melend o et al. 2001). It is the seco nd brightest known 
field RRd star ( Szczygiel & Fabrycky 2007) and it is especially 
suited for high precision CCD photometric observations because 
in its vicinity there are comparison and check stars (GSC 5002- 
0506, GSC 5002-0525, GSC 5002-0566) approximately of the 
same brightness and colour. 

This study presents the first [/-band as well as new BV(RI)c 
time-series data for this RRd star and the fundamental param- 
eters which could be derived from the frequencies. Sect. 2 de- 
scribes the observations and reductions, Sect. 3 presents the fre- 
quency analysis and Fourier fitted colour curves. In Sect. 4 a 
short discussion is given and the conclusions are drawn. 



2. Observations and reduction 

2.1. Observations 

The observational material of V372 Ser was collected in the 
spring of 2007 with two telescopes (see Table Q] for details). 
The 1-m RCC telescope is mounted at Piszkesteto Mountain 
Station of the Konkoly Observatory. It was used with the UV- 
enhancement coating Versarray 1300B camera constructed by 
Princeton Instruments. This device contains a back illuminated 
EEV CCD36-40 1 340 x 1300 chip that corresponds to a 6:6x6^8 
field of view (FOV) with 07303 pixeL 1 resolution. The camera 
was driven by QPAsO, a Tel based observing software running 
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Table 1. Logbook of the observations of V372 Ser. 



HJD-2 400 000 


No. of frames 


quality 9 


Telescope 


54217. 3745-. 5945 


305 


2 


RCC 


54220.4249-.4670 


55 


3 


RCC 


54221.3889-.5685 


250 


3 


RCC 


54222.3767-.5986 


310 


2 


RCC 


54223. 3726-.5865 


300 


2 


RCC 


SZ19A9 S 1 £.1 
JHZ'tZ. JJJJ- .J 1 OZ 




9 


prr 
l\ v v 


54244.3418-.3808 


48 


3 


RCC 


54245. 3444-.5280 


223 


2 


RCC 


54245. 3971-.6290 


211 


1 


IAC80 


54248.3915-.6535 


165 


1 


IAC80 


54249.4283-.6746 


250 


1 


IAC80 


54250.4068-.6728 


270 


1 


IAC80 


54251.4152-.6580 


224 


1 


IAC80 



a Quality of the night is characterized by the extinction coefficient k' v 
in units V mag-(airmass)~' : 1: stable and k' v < 0.25, 2: stable and 
0.25 < k' v < 0.5, 3: k' v > 0.5 or k' v < 0.5 and unstable or k' v > 0.5 
and interruptions because of clouds. 



on Linux platforms developed by Konkoly staff. The IAC8CQ 
telescope of Teide Observatory with its standard camera was the 
other telescope. Its basic parameters are: EEV UV-coated back 
illuminated 2148x2048 chip, 1 1' x 10:5 FOV with 07305 pixel -1 
resolution. Standard Johnson UBV and Kron-Cousins (RI)c fil- 
ters were used in the observations. 

The total number of images is 2812 in all colours, more than 
560 for each filter. The typical exposure times were 150, 80, 40, 
20, 20 s (IAC80), 120-180, 60, 15-40, 5-20, 10-20 s (RCC) for 
U, B, V, Rc, Ic, respectively. Each night sky flats (or alternatively 
dome flats) and some bias frames (5-50) were taken as main cal- 
ibration images. 

2.2. Data reduction & photometry 

We used the iraf/ccdrelQ package for the standard reduction 
procedures: bias and flat field correction. Other corrections 
(e.g. dark current, deferred charge, non-linearity) were also in- 
vestigated, but they were found to be negligible. 

The star brightness was determined by using the aperture 
photometry task daophot/phot of Iraf. This method provides ro- 
bust flux estimates with small errors for isolated stars such as 
V372 Ser, its comparis on and the check stars. For a fi nding chart 
see Fig. 1 in the paper lGarcia-Melendo et all (1200 ll) . Basically, 
we c arried out relative photometr y using GSC 5002-0506 (star 
A of iGarcia-Melendo et alj|200lh as primary comparison star. 
GSC 5002-0525 and GSC 5002-0566 (stars B and C, respec- 
tivel y) were the check stars. The red star GSC 5002-0560 (star 
D of IGarcia-Melendo et alJl200TI) was excluded from our study 
because its colour indices differ from those of the variable and 
the bright ness we measured differs si gnificantly from the value 
given by Garc ia-Melendo et alJ d200lh . It might be a long period 
red variable of small amplitude. 



1 The 0.82m IAC80 Telescope is operated on the island Tenerife by 
the Instituto de Astrofisica de Canarias in the Spanish Observatorio del 
Teide 

2 iraf is distributed by the NOAO, operated by the Association of 
Universities for Research in Astronomy Inc., under contract with the 
NSF. 



2.3. Transformation into the standard system 

The instrumental magnitudes were transformed into the 
UBV(RI)c J ohnson-C o usins system. The standard sequence 
published by iLandolj ( 1 1973b was used as a source of stan- 
dard magnitudes. Our standard UBV(RI)c magnitudes for stars 
A, B, C were determined on the nights HJD-2 454 200=23 
and 48. Within the er r or our values were identical to those of 
IGarcia-Melendo et al.l (|2001|) . therefore, we used the weighted 
mean (giving half weight to our measurements of 2454223). 
Table|2]gives our accepted magnitudes and colour indices. Each 
night colour differences « A - «b, b A - £>b, ■■■ were determined 
from all frames and correlated with Ua - Ub,...,B a - B%,..., 
where u,b, ... are the instrumental magnitudes of stars A, B, C. 
Since almost perfect linear relations and random scatter were 
found the magnitude differences of V372 Ser were transformed 
to the standard U BV(RI)c magnitudes by the formulae 

U. = U A + 1.016(w,-u A ) + 0.012, 
B„ = B A + l-018(K-b A ) + 0.011, 

V, = Va + 1.005(v*-v a ) + 0.001, (1) 
Rc,* = Rc,a + 1.001(r.-r A ) + 0.002, 
h,* = /c,a + 0.981(/«-/ a ) + 0.013, 

where the asterisk indicates V372 Ser. Between 
HJD-2454245=0.3971-.5280 we have simultaneous ob- 
servations with the telescopes RCC and IAC80. Zero point 
corrections -0.043, -.025, -0.034, -0.019, +0.004 were added 
to the RCC U, B, V, (R, I)c magnitudes in order to obtain 
congruence of the RCC and IAC80 light curves. 

The typical errors of individual observations were between 
0.01 - 0.02 mag in B, V, (R, I) c and 0.03 mag in U. 

The construction of our observational databas^l is simple: 
the file includes the heliocentric Julian date (HJD) and magni- 
tudes in the Johnson-Cousins reference U,B,V,Rc,Ic, respec- 
tively. In the last column of the table a quality flag is given, 
it indicates the 13 frames as well which were excluded from 
the analysis because of bad atmospheric conditions (e.g. sudden 
transparency loss exceeding 1.5 mag). 

3. Analysis of the light curves 

3.1. Pulsation frequencies 

We performed standard Fourier frequency analysis on our time 
series data usin g the facilities of the program package MuFrAn 
jKollafh| [T990). The frequencies of the fundamental and first 
overtone modes (f = 2.121841 d" 1 , f\ = 2.851189 d~') and 
their +1,2,... daily aliases were clearly observed. With the first 
approximation of these two frequencies and some harmonics a 
simultaneous nonlinear fit was carried out. Then the light curves 
were whitened out by subtracting the result of the fit. The Fourier 
spectra of the residual light curves contain strong peaks which 
can be identified with the coupling term fo + fi and its aliases. 
Their appearance is an unambiguous signal of strong nonlinear 
coupling in the stellar pulsation. In order to identify further fre- 
quency components, a successive prewhitening procedure was 
performed. The subsequently found frequencies were fitted si- 
multaneously with all previously identified ones and their overall 
contribution was subtracted from the original data set. This iter- 
ative method increases the signal detection probability by elim- 
inating the strong alias structure in the vicinity of a frequency 

3 Available also at the webpage 

http : / /www . konkoly . hu/ staff /benko/pub . html 
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V 


U-B 




B—V 


V-R c 




V-I c 




GSC 5002- 
GSC 5002- 
GSC 5002- 


-0506 
-0525 
-0566 


11. 500 + .024 
13.524 + .026 
12.925 + .025 


0.021 ± 
0.156 ± 
0.052 + 


.005 
.025 
.017 


0.534 + .004 
0.655 ± .005 
0.549 ± .005 


0.336 ± 
0.384 ± 
0.329 ± 


.001 
.004 
.004 


0.664 ± 
0.764 + 
0.668 ± 


.002 
.004 
.008 




50539 



50540 



11 
11.5 



50560 



50569 



54242 



54243 



50541 



50542 



50545 



50548 



50550 



50551 



50552 



50570 



50571 



50573 



50574 



50575 



50576 



50577 



54244 



54245 



54246 



54247 
HJD-2400000 [day] 



54248 



54249 



54250 



50553 



50578 



54251 



50554 









iVr 




50588 



50639 



54224 



54252 



Fig. 1. V light curve of V372 Ser. Blue, green, and red symbols denote observations of iGarcia-Melendo et al. (2001), observations 
with the lm RCC telescope at Konkoly Observatory and the IAC80 telescope at Teide Observatory, respectively. Line: fitted Fourier 
solution Eq. (f2]). 

Table 3. Fourier decomposition of the UBV(RI)c light curves. Units of frequency, amplitude, phase are d -1 , mag, and radians, 
respectively. The first line gives Aqq in different bands. 



Vjk 


frequency 


<> 

jk 


W) 


jk 


(B) 


A™ 

jk 


(V) 


jk 




A ( « 

jk 





fo 

2/o 

3./,, 

./', 

2/i 

3/i 

4/i 

fo+fx 
fi ~ fo 
/o + 2/i 
2/o+/i 
2/o+2/! 
fo + 3/ 
3/o+/ 



2.121840 

4.243680 

6.365520 

2.851188 

5.702376 

8.553564 

11.404752 

4.973028 

0.729348 

7.824216 

7.094868 

9.946056 

10.675404 

9.216708 



11.649 

0.20136 

0.04098 

0.00823 

0.26444 

0.04182 

0.01553 

0.00542 

0.09113 

0.04001 

0.02698 

0.02190 

0.01630 

0.01244 

0.00501 



1.41705 
5.01763 
2.86820 
0.41696 
3.98233 
1.03100 
4.45459 
4.09995 
4.48008 
0.88727 
1.83531 
5.02404 
5.25972 
5.39284 



11.675 

0.19769 

0.03473 

0.00816 

0.26206 

0.04421 

0.01461 

0.00729 

0.08176 

0.05189 

0.02363 

0.02122 

0.01961 

0.01515 

0.01189 



1.33116 
4.88158 
3.48286 
0.44289 
3.99328 
1.27607 
4.52769 
4.16265 
4.56988 
1.01849 
2.01262 
5.27763 
5.21436 
5.45943 



11.264 

0.15399 

0.02771 

0.00828 

0.20591 

0.03652 

0.01220 

0.00632 

0.06282 

0.04279 

0.01894 

0.01600 

0.01453 

0.01196 

0.01086 



1.28978 
4.83961 
3.53757 
0.42275 
3.98314 
1.21576 
4.37625 
4.16672 
4.53610 
1.04737 
1.97154 
5.30904 
5.11548 
5.51522 



11.031 

0.12873 

0.02248 

0.00588 

0.16723 

0.02872 

0.00989 

0.00525 

0.05197 

0.03747 

0.01532 

0.01470 

0.01267 

0.01041 

0.00905 



1.22187 
4.89025 
3.39453 
0.38875 
4.01917 
1.33326 
4.34081 
4.16621 
4.62943 
1.11416 
2.02126 
5.30237 
5.05965 
5.63708 



10.821 

0.09383 

0.01590 

0.00541 

0.12366 

0.02021 

0.00878 

0.00360 

0.03790 

0.02532 

0.01194 

0.01064 

0.00855 

0.00697 

0.00637 



1.15808 
4.81655 
3.80429 
0.32015 
3.93914 
1.46891 
4.21666 
4.14347 
4.56013 
1.08750 
1.96231 
5.30350 
5.12237 
5.68115 



peak, enabling one to identify hidden low amplitude components 
suppressed by aliases and separate frequency peaks that are close 
to each other. As a result, 14 frequencies of the light variation 
were found (see Table[3j. 

To determine the final set of frequencies we investigated 
the CCD V data set obser ved in the years 1997-1998 by 
IGarcia-Melendo et al.l d2001l) . Performing the above described 
frequency search we obtained a Fourier solution with 12 ele- 
ments of very similar parameters to our previous ones, both for 
the frequencies (f = 2.12250 d" 1 , fx = 2.850833 d~') and the 



amplitudes. Differences between periods are 3.78 and 12.64 sec- 
onds for /j" 1 and , respectively: the frequency content of the 
light variations seems to be unchanged within the last ten years. 
To verify this conjecture we tried to make a common fit of all 
the available V data. We remark that because of the lack of data 
for 10 years to unify and handle the two data sets in a common 
Fourier analysis is not practical due to the very unfavourable 
window function. None of the further determined Fourier param- 
eters (frequencies, amplitudes and phase) based only either of 
data sets gave us an acceptable fit for the unified data. Therefore, 
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we refined the frequency analysis by a successive approximation 
on the combined data, where the residuals of the non-linear fit 
were minimized while the structure of the decomposition (num- 
ber of harmonics and linear combinations) was fixed and the two 
frequencies were independently changed in their limiting inter- 
vals. As a result of this process we could fit the total light curve 
at the frequencies f = 2.121840 d _1 , /i = 2.851188 d l with 
the rms accuracy of 0.014 mag. 

Since the frequencies were found to be stable we can em- 
ploy the spectrum aver aging method (SAM) as implemented by 
iNagv & Kovacsl (120061) to find additional peaks which remained 
hidden in the individual spectra. The two spectra were added 
weighted inversely with their variances. The morphology of the 
peaks of the summed spectrum confirmed again that the fre- 
quencies were really constant within their errors. We subtracted 
a non-linear fit of the accurate frequencies and their harmon- 
ics from both light curves separately, then the Fourier spectra of 
residuals were summed again to find the next frequency and so 
on. At the end of this process we achieved a spectrum without 
any significant additional peaks at 6cr level (for the definition 
of significance level of a peak see lAlcock et al.1 12000). Our con- 
clusion is that additional signal has not been found: the Fourier 
decomposition of our data revealed all significant frequencies. 

3.2. Fourier amplitudes and phase 

Using the obtained frequencies, their harmonics and linear com- 
binations non-linear fits were computed determining amplitudes 
and phase for all bands. We assumed the light curves m(f) in the 
form 

m(t) = A$ + J] A jk ic) sm[2nv jk (t - to) + <pf], (2) 

jk 

where c = U,B,V,R c Jc, v jk = jfo + kf u j;k = 0,+l,±2, 
and the epoch of an arbitrarily chosen phase was ?o = 2450539. 
The results are given in Table [3] These parameters allowed us to 
construct synthetic light curves. As an example the synthetic and 
observed V light curves are shown in Fig. Q] The residual rms 
scatters of the fits are 0.012, 0.015, 0.015, 0.018 for B, V,R c Jc 
bands, respectively, and it is 0.028 for U. These values of scatter 
are compa tible with the observational ac curacy. For the unified V 
filter data dGarcia-Melendo et af. 2001 and ours) we also found 
a fit with rms=0.014, similar amplitudes, and phase as given in 
Table [3] It indicates, that on 10 years scale not only frequencies 
and amplitudes were stable but phase as well. Folded V light 
curves are given in Fig [2] 

The amplitudes of all detected frequencies increase towards 
shorter wavelengths and the amplitudes of the linear combina- 
tion frequency components are smaller than those of their con- 
stituent frequencies. All of the detected coupling terms are at 
positive linear combination frequencies of the modes except for 
fo - f\ . This seems t o be in agreement with model results of 
lAntonello & Aikawal d 1998b . According to their non-linear tran- 
sient double-mode Cepheid models, frequency terms of positive 
linear combinations have always larger amplitudes than those of 
negative linear combination. 

To estimate the accuracy of the parameters in Table [3] Monte 
Carlo simulations were carried out. 1000 artificial time series 
were constructed for data in each band with the help of the de- 
termined Fourier parameters, Gaussian noise and sampling ac- 
cording to the real observations. The epochs of the beginning of 
data sampling were randomly chosen for the first trial then they 
were shifted by 0.01 day for each subs equently produced syn- 
thetic time series. As it was stressed bv lBenko et al.l d2006l) this 



10.5 




phase 

Fig. 2. Folded V light curve of V372 Ser, blue and red symbols 
denote observations of Garcia-Me lendo et alj d2001l) and the ob- 
servations with the lm RCC and IAC80 telescopes, respectively. 
Upper panel: folded with Pi, middle panel: prewhitened with 
Po, folded with P\, lower panel: prewhitened with P\, folded 
with Pq. 



is important for the interpretation of possible amplitude varia- 
tions such as mode changes. The accuracies were found to be 
between 0.0025 and 0.0007 mag for all amplitudes. The highest 
errors are for the filter U data and the smallest ones are for fil- 
ters Rq and Ic- Typical errors of the phase are between 0.01 and 
0.004 radians for the two main frequencies. 

4. Discussion and conclusions 

Table [4] summarizes the relevant data of V372 Ser derived 
from our Fourier analysis. The average of the amplitude ratios 

A o"i7 A io = 1-319 + .008, m = U,B, V,R c Jc shows that the am- 
plitude ratio is independent of colours at 3<x level. 

From the point of view of frequencies and amplitudes 
V372 Ser is a double mode RR Lyrae star: the period ratio 
0.7441951 is in the canonical range give n by non-linear theo- 
retical values expected for RR Lyrae stars dCox et all l983). The 
first overtone is the dominant mode. In different bands identical 
frequency content has been found: the atmosphere is grey from 
the point of view of transmitting pulsation frequencies. 

The magnitude-averaged colour indices obtained by fitting 
the date with Eq. d2]l are given in Table |H they differ slightly 
from those of iGarcia^Melendo et al.l d200 Ik (B-V) = 0.38, (V- 
R c ) = 0.26, except for (V - / c > = 0.57. We attribute this differ- 
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Table 4. Summary of the relevant data of V372 Ser. 



/>,(days) 




0.3507310 + .0000001 


Po(days) 




0.4712891 ± .0000002 


Pi/Po 




0.7441950 






11.264 


(U-B) = A$- 


- A m 


-0.026 


(B-V)=A®- 


^00 


0.411 


(V-R C )=A<£ 


- A (Rc ' 
A oa 


0.233 


(V-Ic)=A™ 


- A i c) 


0.443 


A m 




0.2059 


A w 




0.1534 


A W) , a (U) 
A ()i I A 10 




1.313 


a (b), a (B) 
A ni > A ip 




1.326 


A (v) /A m 

A 0l ' A 1Q 




1.337 


A (Rc, /A (Rc> 

A 0) ' A 10 




1.299 


A 01 ' A 10 




1.318 



u 




U 




U 




11 




11 


- ^ - 


11 




11.2 


- f\ - 


11.2 




11.2 




11.4 




11.4 




11.4 




11.6 
11.8 




11.6 
11.8 




11.6 
11.8 


: \ ! ] 


12 


- w - 

" " 


12 


1 _ 
' " 


12 


' " 




0.4 0.6 


8 


0.4 0.6 


8 


0.4 0.6 




HJD-2454249 




HJD-2454250 




HJD-2454251 



Fig. 3. Observed U light curves of a medium, large, and small 
ascending branch 



ence to their much smaller sampling of the B,Rq,Ic bands and 
to a possible misprint in their Rq - Ic = 0.29 instead of 0.19. 
V372 Ser has reddish colours even if they are dereddened by 
the upper limit E(B - V) = .08 for the galactic coordinates of 
V372 Ser (ISchleeel et al.lfl998l) . 

Three U light curves of night quality 1 are reproduced in 
Fig. [3] They are smooth: hump or bump indicating atmospheric 
shocks are not present, in contrast with the ascending b ranch 
of e.g. SU Dra jPreston & Paczvnski|[l964l;lBarczal l2002). The 
duration of an ascending branch is a Pi/2. 

Upper limits have been estimated for secular change P = 
dP/dt of the periods from Monte Carlo simulations and a change 
10~ 6 in the frequencies in order to find minimal standard devia- 
tion of the V synthetic light curve and observations: 

|Pol < 6.0 X 10" 1 'day/day, 

\Pi\ < 2.4 x 10" 1 'day/day. 
For comparison we mention that |P| < 7 x 10~ n is typical during 
the whole lifetime of an RRL star pulsating in single mo de ex- 
cept f or the short phase before leaving the instability strip ( Smith 
1 1995b . Upper limits for secular amplitude and phase changes are 
OidAv/dt) ~ 10- 6 mag/day, 0(d<p Y /dt) ~ lO^radians/day. 

The ratio Pi/P = 0.744195 and P = 0.471289 are ap- 
proximately equal to those of the RRd stars of M3 and IC 
4499, both are moderately metal deficient globular clusters 
([M] = -1.57,-1.50, r espectively). Interp olating Pi/P ,P in 
the Petersen diagram of [Bono et alj {1996) (derived from non- 
linear, non-local, time dependent models with the new opacities) 



gives M — 0.65vVfo and luminosity logL a 1.72. Turbulent 
convection and more sophisticated hydrodynamic treatment 
are a dded to the above model assumptions in dSzabo et al.l 
2004). Interpolation in their diagrams and tables gives [M] » 
-IJ,M >O.7M and logLa 1.60. 
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